arXiv:1503.05383vl [math.PR] 18 Mar 2015 


Modem Stochastics: Theory and Applications 1 (2014) 167-180 
DOI: 10.15559/15-VMSTA18 


Practical approaches to the estimation of the ruin 
probability in a risk model with additional funds 


Yuliya Mishura^, Olena Ragulina^*, Oleksandr Stroyev'’ 

‘^Taras Shevchenko National University of Kyiv, 

Department of Probability Theory, Statistics and Actuarial Mathematics, 

64 Volodymyrska Str., 01601 Kyiv, Ukraine 
'°Yuriy Fedkovych Chernivtsi National University, 

Department of Mathematical Modelling, 

2 Kotsjubynskyi Str, 58012 Chernivtsi, Ukraine 

myus@univ.kiev.ua (Yu. Mishura), ragulina.olena@gmail.com (O. Ragulina), 
o.stroiev@chnu.edu.ua (O. Stroyev) 

Received: 4 January 2015, Accepted: 19 January 2015, Published online: 2 February 2015 

Abstract We deal with a generalization of the classical risk model when an insurance com¬ 
pany gets additional funds whenever a claim arrives and consider some practical approaches 
to the estimation of the ruin probability. In particular, we get an upper exponential bound and 
construct an analogue to the De Vylder approximation for the ruin probability. We compare 
results of these approaches with statistical estimates obtained by the Monte Carlo method for 
selected distributions of claim sizes and additional funds. 
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1 Introduction 

Let (17,5^, P) be a probability space satisfying the usual conditions, and let all the 
objects be defined on it. We deal with the risk model that generalizes the classical one 
and was considered in [10]. 
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In the classical risk model (see, e.g., [1, 7, 11]), an insurance company has an ini¬ 
tial surplus X > 0 and receives premiums with constant intensity c > 0. Claim sizes 
form a sequence (fi)i>i of nonnegative i.i.d. random variables with c.d.f. Fi{y) — 
^ y] finite expectation = pi. The number of claims on the time 

interval [0, t] is a homogeneous Poisson process {Nt)t>o with intensity A > 0. 

In addition to the classical risk model, we suppose that the insurance company 
gets additional funds rji when the jth claim arrives. These funds can be considered, 
for instance, as additional investment income, which does not depend on the surplus 
of the company. We assume that (r 7 i)i>i is a sequence of nonnegative i.i.d. random 
variables with c.d.f. F 2 {y) = P[r 7 i < y] and finite expectation £[ 771 ] = p 2 - The 
sequences {r]i)i>i and the process {Nt)t>o are mutually independent. Let 

(5'i)t>o be the filtration generated by (77i)i>i, and {Nt)t>o. 

Let Xt{x) be the surplus of the insurance company at time t, provided that its 
initial surplus is x. Then the surplus process {Xt{x))t>o is defined as 

Nt 

Xt{x) = X + ct- r]i), t>0. (1) 

i=l 

Note that we set X]i=i(Ci ~ Vi) = ^ in (1) if = 0. 

The ruin time is defined as 

t{x) = inf{f > 0: Xt{x) < O}. 

We suppose that t(x) = 00 if Xt{x) > 0 for all f > 0. The infinite-horizon ruin 
probability is given by 

tpix) = P[inft>o Xtix) < 0 ], 

which is equivalent to 

ipix) = P[r(x) < cx)]. 

The corresponding infinite-horizon survival probability equals 

'•p{x) = 1 — ^^(x). 

Note that the ruin never occurs if P[^i — Pi < 0] = 1. If P[^i — Pi > 0] = l,then 
we deal with the classical risk model. So in what follows, we assume thatP[^j — Pi > 
0] > 0 and P[^i — p^ < 0] > 0. In this case, if c — \y,i + X 112 < 0, then ip{x) = 0 for 
all a; > 0; if c — Xfii + Xy ,2 > 0, then lima;_>,_|_oo ‘fix) = 1 (see [10, Lemma 2.1]). 

In this paper, we consider some practical approaches to the estimation of the ruin 
probability. In particular, we get an upper exponential bound and construct an ana¬ 
logue to the De Vylder approximation for the ruin probability. Moreover, we compare 
results of these approaches with statistical estimates obtained by the Monte Carlo 
method for selected distributions of claim sizes and additional funds. 

Paper [10], where this risk model is considered, is devoted to the investigation of 
continuity and differentiability of the infinite-horizon survival probability and deriva¬ 
tion of an integro-differential equation for this function. When claim sizes and addi¬ 
tional funds are exponentially distributed, a closed-form solution to this equation can 
be found. 
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Theorem 1 ([10], Theorem 4.1). Let the surplus process {Xt{x))t>o follow (1) under 
the above assumptions, the random variables and rji, i > 1, be exponentially 
distributed with means pi and p 2 correspondingly, and c — Xpi + Xp 2 > 0. Then 


<p{x) = 1 + 7 — 
[ca 


Xpi{l - ap2) _ 

A)(l - ap2){pi + P 2 ) + Xp2 


( 2 ) 


for all X > 0, where 


_ XpiP2 + Cpi — CP2 — y/c^ifJ-l + T 2 ) + + ‘^CpiP2{c — Xpi + Xp2) 

2cpiP2 

Remark 1 ([10], Remark 4.1). It is justified in the proof of Theorem 1 that a < 0 
and 

^ _ Xpi{l-ap2) _ ^ Q 

{ca — A)(l — ap 2 ){Ti + M 2 ) + Xp 2 

So the function p{x) defined by (2) satisfies all the natural properties of the survival 
probability. In particular, this function is nondecreasing and bounded by 0 from below 
and by 1 from above. 

It is well known that even for the classical risk model, there are only a few cases 
where an analytic expression for the survival probability can be found. So numerous 
approximations have been considered and investigated for the classical risk model 
(see, e.g., [1, 2, 4, 5, 7, 11]). “Simple approximations” form a special class of approx¬ 
imations for the ruin or survival probabilities. They use only some moments of the 
distribution of claim sizes and do not take into account the detailed tail behavior of 
that distribution. Such approximations may be based on limit theorems or on heuristic 
arguments. The most successful “simple approximation” is certainly the De Vylder 
approximation [5], which is based on the heuristic idea to replace the risk process 
with a risk process with exponentially distributed claim sizes such that the first three 
moments coincide (see also [7, 11]). This approximation is known to work extremely 
well for some distributions of claim sizes. Later, Grandell analyzed the De Vylder 
approximation and other “simple approximations” from a more mathematical point 
of view and gave a possible explanation why the De Vylder approximation is so good 
(see [8]). 

We deal with the case where the claim sizes have a light-tailed distribution. The 
rest of the paper is organized as follows. In Section 2, we get an upper exponential 
bound for the ruin probability, which is an analogue of the famous Lundberg inequal¬ 
ity. Section 3 is devoted to the construction of an analogue of the De Vylder approxi¬ 
mation. In Section 4, we give a simple formula that relates the accuracy and reliability 
of the approximation of the ruin probability by its statistical estimate obtained by the 
Monte Carlo method. In Section 5, we compare the results of these approaches for 
some distributions of claim sizes and additional funds. Section 6 concludes the paper. 


2 Exponential bound 

To get an upper exponential bound for the ruin probability, we use the martingale 
approach introduced by Gerber [6] (see also [3, 7, 11]). 
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Let 


Nt 


Ut=ct- - Vi), t > 0- 


2=1 


For all i? > 0, we define the exponential process (Vt (i?))t>o by 

Vt{R)=e-’^^\ 

Lemma 1. If there is R > 0 such that 


^ + oo 


^+oo 


A e^y dFi (y) • dF 2 (y) - 1 j = cR, 

then {Vt{R))t>o is an {'^t)-tntirtingale. 

Proof. For all i? > 0 such that E[e^^‘] < oo, if any, we have 


(3) 


E[yt(i?)] =e“^'^*E 


Nt 


exp< R'Y^i^i - Vi) 


cm 


E 

j=o 


J'- 




(4) 


= exp{f(AE[e^^^‘ — A — ci?)}. 

If there is i? > 0 such that (3) holds, then E[e^^’] < oo, and for all t 2 > ti > 0, 
we have 


E[Vt,{R)/SH] =E 


Nt. 


= E 


exp' 


exp' 


ct2 - - Vi) 

\ i=l / 

/ ^*1 N 

Cti - J2(^i - V^) 




i=l 


X E 


Nt. 


exp-^ -A I c{t2 -ti)- {Ci - Vi) 

I V i=Nt^ / 


S'*! 


= E[Ft,(i?)] -E 

= E[F*,(i?)]. 

Here we used the fact that 


Nt 


*2-*! 


exp<^c(f2 - ft) - i^t-Vt) 


E 


Nt. 


exp|-^|^c(f2 - ft) - 

= exp{(f2 — fi)(AE[e^^^‘“’’*^] — A — c-R)} = 1 


by (3) and (4). 

Thus, {Vt{R))t>o is an -martingale, which is the desired conclusion. 


□ 
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Theorem 2. If there is R > 0 such that (3) holds, then for all x >Q, we have 

f’ix) < e-^^. (5) 

Proof. It is easily seen that t{x) is an (5^t)-stopping time. Hence, t{x) A T is a 
bounded (5^t)-stopping time for any fixed T > 0. The process (Vt(J?))t>o is an 
(5^t)-martingale by Lemma 1. Moreover, {Vt{R))t>o is positive a.s. by its definition. 
Consequently, applying the optional stopping theorem yields 


= ]E[ 14 ( 2 ,)(.R) • I{t(x)<t}] -f ]E[Lr(.R) • I{r(x)>T}] 
> E[14(2,)(.R) ■ I{t(x)<t}] 


N, 


(x) 


= E exp|-.R I^ct(x) - ^ (Ci - Vi) 
>e^“-P[T(a:) < T] , 

where is the indicator of an event. This gives 

F[t{x) <T]< 

for all T > 0. Letting T —oo in (6) yields 

P[r(x) < oo] < e~^^, 


■I 


{r(x)<T} 


(6) 


which is our assertion. □ 

Example 1. Let the random variables and p,, z > 1, be exponentially distributed 
with means pi and /i 2 , respectively. Then (3) can be rewritten as 

Af"-^ 

V(l - pii?)(l + P2-R) 
where R G (0, l/p-i). This condition is equivalent to 

cpi/r2.R^ + (Apip2 + c/zi — cfi2)R^ — {c — A/zi -I- \fL2)R = 0. (7) 

If c — A/ri -f Ap 2 > 0, then 

+ nl) + + 2cpip2(c - A^i -I- XfL2) > 0. 

So there are three real solutions to (7). They are 

Ri = 0 , 

- _ XfJ.ltl2 + Cfj.1 - Cfj.2 - s/A{c, X, Hi, fj,2) 

' “ 2cpiP2 ’ 
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where 


Rs 


+ C/Ti - C^X2 + v^A(c, A,/ri,/X2) 

2c^1^2 


A{c,X,fj.i,fi2) = + M2) + ^^MiM2 + 2 c/iiM 2 (c - A/ii + Am2)- 

Furthermore, it is easy to check that, in this case, 

|AmiM 2 + CMi - c/i2| < \/ 24 (c, A,mi,M2)- 

From this we conclude that i ?2 > 0 and R 3 < 0. Since 

^(c, A,Ml,M2) < (AM1M2 + CM1 +CM2)^, 


we have 

n (AmiM 2 + CMi + CM 2 ) - (AmiM 2 + CMi - CM 2 ) 1 

R 2 < -^- < —• 

2cmiM 2 Ml 

Hence, R 2 is a unique positive solution to (7), and an exponential bound (5) can be 
rewritten as follows: 

^{x) < (8) 

Comparing (8) with (2), we see that the exponential bound and the analytic expression 
for the ruin probability differ in a constant multiplier only. 

If c — Ami + Am 2 < 0, then > ^ 2 , which gives AmiM 2 + CMi — CM 2 > 0. 
Let R 2 and R 3 be two nonzero solutions to (7). Applying Vieta’s formulas yields 
R 2 + R 3 < 0 and R 2 R 3 > 0. Consequently, if R 2 and R 3 are real, they are negative. 
Thus, (7) has no positive solution, and Theorem 2 does not give us an exponential 
bound for the ruin probability. Indeed, in this case, tpix) = 1 for all a: > 0 (see [10, 
Lemma 2.1]). 


3 Analogue to the De Vylder approximation 


To construct an analogue to the De Vylder approximation, we replace the process 
{Ut)t>o with a process {Ut)t>o with exponentially distributed claim sizes such that 

£[(//=] =E [(//=], fc = l,2,3. (9) 

Since the process {Ut)t>o in this risk model is determined by the four parameters 
(c, A, Ml, M 2 ) in contrast to the classical risk model, where it is determined by three 
parameters, we use the additional condition 


m m 

M2 M2' 


( 10 ) 


Note that we could have used the condition = E[{7)^] instead of (10), but it 

would have led to tedious calculations and solving polynomial equations of higher 
degree. 
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Let be a sequence of i.i.d. random variables exponentially distributed with 

mean /xi. Similarly, let be a sequence of i.i.d. random variables exponentially 

distributed with mean /x2. An easy computation shows that 


and E[f,f] = 

Let E[^f] < oo and E[77f] < oo. Then we have 

■ Nt 


( 11 ) 


E[17t] = cf — E 




2=1 


= ct- AfE[^i - 



' Nt 


r / Nt \ 2-| 

E[Uf] = (ctf-2ctE 

.2=1 

-fE 

- 1 

1 

_1 


= (cf)^ - 2ct ■ AtE[^i - rji] + AfE[(^^ - 77^)^] + (At)^(E[^^ - 77^])^ 
= AfE[(^j - 77i)^] + (ct - AfE[^j - 77i])^ 

= AtE[(C.-770'] + (E[^7t])^ 


E[Uf] = {ctf -3{ctfE 


Nt 




,i=l 


3ctE 


' Nt 




-E 


Nt 


3n 




= {ctf - 3 {ctf' ■ AfE[^i - 77*] + 3 cf(AfE[(^i - 77^)^] 

+ {Xtf{E[^t - 77.])" - AtE[(e. - 77,)"] - 3 (At) 2 E[(^, - 77 ,)'] 

+ AfE[(^, - r,tf] E[(^, - Tjt)2 - (Af)"(E[e. - 77 ,])' 

= -AfE[(^j - T]t )2 + {ct - AfE[^i - T]t]f 
+ 3 AfE[(^i - 77 i)^] {ct - AfE[^i - rjt]) 

= -AtE[(e. - v^f] + (E[Lt])' + 3(E[172] _ (E[C/t])') E[Ut]. 


Applying similar arguments to the process {Ut)t>o, we conclude that (9) is equiv¬ 
alent to 

{ ct - AfE[^^ - 77i] = ct- AfE[|i - fji], 

XtE[i^t-V^f] =AfE[(|,-77,)2], (12) 

-AfE[(e.-p.)"] =-AfE[(f,- 77 ,) 3 ], 

By (11) we can rewrite (12) as 


{ c - A(77 i - 7x2 ) = c - A(/2i - 7x2), 

AE[(^j - 77z)^] = 2A(/i? - A 1 /I 2 + A 2 ): 

AE[(^j - Vif] = 6A(7i? - + itipi - P-D- 


( 13 ) 
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Substituting /i 2 = /i2Mi/Mi into the second and third equations of system (13) 
yields 


AE[(e.-770'] =2A/l?fl-^ + 4Y (14) 

V Ml Ml/ 

AE[(C. - mf] = 6A/i? (l-^ + 4-4)- (15) 

V Ml MI Ml/ 


Dividing (15) by (14) gives 


Ml (Mi - MiM 2 + Mi) - mf] 

3(Mi - M?M 2 + MimI - M 2 ) ’ 


(16) 


Consequently, we have 

~ ^ M2(Mf - M1M2 + mI) lE[fe - mf] 

3(/rf - + MimI - M 2 ) - Vi)'^]' 

Substituting (16) into (14), we get 


(17) 


, 9A(Mf - + MiMi - Mi)^ (E[(6 - 

2(m? - M 1 M 2 + mD^ (E[(C» - M*)^])^ 


(18) 


Substituting (16)-(18) into the first equation of system (13), we obtain 
3 (mi - M?M2 + MimI - mI) (lE[te - t?i)^])^ 


= c - A(/ri - M 2 ) 1 - 


2(m? - M 1 M 2 + mD^ ]E[(Cz - 


(19) 


Note that since Fi (u) and Foiv) are known, it is easy to find E[(C- — and 

E[(C.-7.)^]if]E[Cf]<-andE["3]<oo. 

By (16) and (17), fli and fl 2 are positive, provided that 


(m? - M?M 2 + M 1 M 2 - M 2 ) E[(?i - Vif] > 0. (20) 


If (20) holds, then pi ^ M 2 - So A is also positive. Moreover, c is positive, provided 
that 


c - A(mi 



3(Mf - MfM2 + MiMi - Mi) (E[(g. - i?.)^])^ A ^ „ ,,,, 

2(M?-MiM2 + Mi)^E[(^,-p,)'] ■ 1 1 


Thus, we get the following result. 

Proposition 1 (An analogue to the De Vylder approximation). Let the surplus pro¬ 
cess (Xt(x))t>o follow (1) under the above assumptions, c — Xp-i + A/i 2 > 0, 
E[^3] 00 ^ E[r 7 f] < c», and let conditions (20) and (21) hold. Then the ruin proba¬ 

bility is approximately equal to 


fjovix) 


_ Xpiiap.2 - 1)_ 

{ca - A)(l - a/i2)(Mi + M 2 ) + Xfl 2 


for all X > 0, where 
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A/ii/22 + cjli — cfl2 — J + Ai) + A^AiAi + 2cAiA2(c — AAi + AA 2 ) 
a = - - - _ - 

2c/Zi/i2 

and the parameters (c, A, Aij A 2 ) ore defined by (16)-(19). 

Remark 2. Note that a < 0 and 

_ AAi(aA2 - 1) _^ Q 

(ca - A)(l - ap.2){fii + A 2 ) + AA 2 

in Proposition 1. Indeed, the parameters (c, A, pi, A 2 ) positive. Moreover, since 
c — Api + Xp 2 > Oj we have c — Api + AA 2 > 0 by the first equation of system (13). 
Hence, Theorem 1 and Remark 1 give us the desired conclusion. 

Remark 3. If claim sizes and additional funds are exponentially distributed, then it 
is easily seen from (16)-(19) that ip{x) = ipDv{x). 

4 Statistical estimate obtained by tbe Monte Carlo method 

Let N be the total number of simulations of the surplus process Xt{x), and let ip{x) 
be the corresponding statistical estimate obtained by the Monte Carlo method. To get 
it, we divide the number of simulations that lead to the ruin by the total number of 
simulations. 

Proposition 2. Let the surplus process {Xt{x))t>o follow (1) under the above as¬ 
sumptions. Then for any e > 0, we have 

P[|V'(a;) — y{x)\ > e] < ( 22 ) 

The assertion of Proposition 2 follows immediately from Hoeffding’s inequality 
(see [9]). 

Remark 4. Formula (22) relates the accuracy and reliability of the approximation of 
the ruin probability by its statistical estimate obtained by the Monte Carlo method. 
It enables us to find the number of simulations N, which is necessary in order to 
calculate the ruin probability with the required accuracy and reliability. An obvious 
shortcoming of the Monte Carlo method is a too large number of simulations N. 
In all examples in Section 5, we assume that e = 0.001 and ^ = 0.001. 

Consequently, N = 3 800 452. 


5 Comparison of results 

5.7 Erlang distributions for claim sizes and additional funds 
Let the probability density functions of and rji be 

= -!- and = f- 

tA^{ki-l)\ lA^{k2-l)\ 

for ?/ > 0, respectively, where ki and k 2 are positive integers. 
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In what follows, hi{R) and h 2 {R), where R> 0, denote the moment generating 
functions of and rji, respectively, that is, 

hi{R) = and h 2 {R) = E[e^’^']. 

An easy computation shows that 

r+oo / ^ 


h2{R) = dF2{y) = 

Moreover, for all i? > 0, we have 

r*+oo 


ki — fiiR 

k2 
k2 - 


0 < i? < —, 
Ml 


fe2 


0 < i? < 


M2 


^ e-^y dF 2 {y) = 


M2-R 


^2 


Thus, condition (3) can be rewritten as 
A ' " ' 


^2 


— A “h C^, 

, k\ — (j,^R J V ^2 3“ y 

where 0 < R < kijfii. Furthermore, we have 

= ^l(O) = Ml. Hdi] = ^ 2 ( 0 ) = M2, 


(23) 


= h'liO) = 


(fci + l)yl 

ki 


{k2 + l)yl 


E[Cf] = + 


kl 


E[ri]=h”{{))= , 

E[,3]=M"'(0)=(^i±lK^±^. 

^2 


Hence, we get 

E[(C. - drf] = Efc'] - 2E[e.]E[p,] + E[p2] 


(fci + l)yl 


— 2/il^2 + 


{k2 + 1)M2 


kl k2 

E[(e. - = E[Cf] - 3E[C"]E[,7,] + 3E[e.]E[p2] _ E[ry3] 

{kl + l)(fci + 2)/if 3(A:i + \)ii\y,2 


kl 


kl 


3(^2 + 1)M2M1 (^2 + 1)(^2 + 2)/i2 


k2 

Substituting E[(^i — and E[(^i — rjiY] into (16)-(19), we obtain the parameters 
(c,A,/ii,/i 2 ). 

Example 2. Let c = 10, A = 4, /ii = 2, /12 = 0.5, fci = 3, k 2 — 2. Then 
R « 0.349093, which may not be an unique positive solution to (23), and 

iPdv{x) = 0.612268 

The results of computations are given in Table 1 . 
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Table 1. Results of computations: Erlang distributions for claim sizes and additional funds 


X 

fix) 

ipDvitr) 


) ■ 100% 

g-fe 


) ■ 100% 


^ Aix') 

0 

0.634149 

0.612268 

-3.45% 


1.000000 

57.69% 


1 

0.492768 

0.439087 

-10.89% 


0.705327 

43.14% 


2 

0.355769 

0.314891 

-11.49% 


0.497487 

39.83% 


5 

0.137737 

0.116142 

-15.68% 


0.174564 

26.74% 


10 

0.023224 

0.022031 

-5.14% 


0.030473 

31.21% 



5.2 Hyperexponential distributions for claim sizes and additional funds 
Let 


^i(y) = Pi.iFiAv) + Pi,2Fi.2{y) H-y > 0, 


where ki > 1, pij > 0, Pip — 1’ Sjli Pi,j di,j = Fij is the c.d.f. 

of the exponential distribution with mean yi j l 

F2iy) = P2,lF2p{y) + P2,2F2Ay) H-hP2.fe2-p2,fe2(y)) v > 0 , 

where k 2 > 1, P 2 ,j > 0, P^j = 1, J2jLi P 2 ,j P 2 ,j = P 2 , and F 2 j is the c.d.f. 
of the exponential distribution with mean fL 2 ,j- 


It is easy to check that 

fei 


hi{R)=Yl 


Pi, 


^ 1 - pijR' 


0 < R < min 


I—— ^ 

LMi,i Ml ,2 Ml,All 


k2 


h2{R) = Y. 


P2d 


iW 


Furthermore, for all i? > 0, we have 


0 < i? < min 




1 1 


1 


{P2,l p2,2 P2,k2 


f+oo h2 

' ^-Ry -Sr^ P'^d 


'dF2{y)=Y 


Hence, condition (3) can be rewritten as 


fc2 




P2j 


1 + /i2ji? 


xiy^h^.y 

\~A-pi,iR ~A + P2 ,jR, 


— A -|- cR. 


where 0 < R < l///i^ 2 , • ■ •, 1/Mi,fci}- Moreover, we have 


k2 


(24) 


E[^i\ = h[{0) = Ypij Pl,j = pi, EA] = AiO) = yP2,j P2,j = P2, 


J=1 




E[er] = Kio) = y 2 pi., py e[A] = h'm = y 2^2,, py 

i=i i=i 
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Table 2. Results of computations: hyperexponential distributions for claim sizes and additional 
funds 


X 

fix) 

1pDv{x) 

flpDv(^) ^ 

) ■ 100% 

g-fe 


) ■ 100% 


' litix) 

0 

0.647560 

0.581428 

-10.21% 


1.000000 

54.43% 


1 

0.540924 

0.521454 

-3.60% 


0.895291 

65.51% 


2 

0.488597 

0.467667 

-4.28% 


0.801545 

64.05% 


5 

0.346390 

0.337363 

-2.61% 


0.575201 

66.06% 


10 

0.202323 

0.195749 

-3.25% 


0.330856 

63.53% 


20 

0.067802 

0.065903 

-2.80% 


0.109466 

61.45% 


25 

0.038194 

0.038239 

0.12% 


0.062965 

64.86% 



fci k2 

IE[C?] = hno) = Kv!] = K'iO) = E 

7 = 1 7 = 1 

Consequently, we get 

/ fcl fc2 \ 

E[(C* - Vif] = 2 ( Epi.7 /^1.7 “ ^1/^2 + E^2j /72j ) . 

\ 7=1 7=1 / 


C fei fcl fc2 fc2 \ 

Ep1'7 ^IJ “ ^2 Ep 1.7 + Ml mIj - E^’2.7 ) • 

7 = 1 7 = 1 7 = 1 7 = 1 / 

Example 3. Let c = 10, A = 4, /7i = 2, ^2 = 0.5, ki = 3, fc 2 = 2, pi_i = 0.4, 
Mi,i = 0.5, pi_2 = 0.3, fii^2 = 2, = 0.3, = 4, p2,i = 0.75, ^ 2,1 = 0.4, 

4 * 2,2 = 0.25, /J, 2,2 = 0.8. Then i? Ri 0.110607, which may not be a unique positive 
solution to (24), and 

ijjDvix) = 0.581428 

The results of computations are given in Table 2. 

5.3 Exponential distribution for claim sizes and degenerate distribution for addi¬ 
tional funds 

Let be exponentially distributed with mean pi and E[? 7 i = P 2 ] = 1- Then condi¬ 
tion (3) can be rewritten as 

-:r — A -f c77, 

1-piR 

where R G (0, 1/pi), which is equivalent to 

+ (c-A/ri)i? +A. (25) 

If c — A/ii + Xp 2 > 0, then it is easy to check that (25) has a unique solution 
R G (0,1/ pi). 
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Table 3. Results of computations: exponential distribution for claim sizes and degenerate dis¬ 
tribution additional funds 


X 

■fix) 



) ■ 100% 

g-fe 


) ■ 100% 


'' ihix') 

0 

0.637998 

0.582498 

-8.70% 


1.000000 

56.74% 


1 

0.549737 

0.482780 

-12.18% 


0.822610 

49.64% 


2 

0.465171 

0.400133 

-13.98% 


0.676687 

45.47% 


5 

0.277026 

0.227808 

-17.77% 


0.376678 

37.97% 


10 

0.113399 

0.089093 

-21.43% 


0.141886 

25.12% 



Since 

E[e.]=/71, E[C,']=2^?, E[c3]^g^3^ 

E[?7i] = ^ 2 , E[? 72 ] = ^2^ E[r]f] = fil, 

we get 

E[(6 - V^f] = 

E[(Ci - Vi)^] = ^lA - 6MiM 2 + 3^iAi2 - lA- 
Example 4. Let c = 10, A = 4, /ii = 2, /j ,2 = 0.5. Then R ^ 0.195273 and 
ij;Dy{x) = 0.582498 

The results of computations are given in Table 3. 

6 Conclusion 

Tables 1-3 provide results of computations when the initial surplus is not too large. 
In this case, the statistical estimates obtained by the Monte Carlo method can be 
used instead of the exact ruin probabilities to compare an accuracy of the exponential 
bound and the analogue to the De Vylder approximation. To get appropriate statis¬ 
tical estimates by the Monte Carlo method for large initial surpluses, the number of 
simulations must be exceeding. The results of computations show that the exponen¬ 
tial bound is very rough. The analogue to the De Vylder approximation gives much 
more accurate estimations, especially in the case of hyperexponential distributions 
for claim sizes and additional funds. Nevertheless, it is heuristic, and its real accuracy 
is unknown. 
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